Army  Research  Laboratory 


111* 


A  One-Dimensional 
Atmospheric  Boundary  Layer 
Model:  Intermittent  Wind 
Shears  and  Thermal  Stability  at 

Night 


Arnold  Tunick 


ARL-MR-494 


June  2001 


Approved  for  public  release;  distribution  unlimited. 


20010718  078 


The  findings  in  this  report  are  not  to  be  construed  as  an 


designated  by  other  authorized  documents. 


Citation  of  manufacturer's  or  trade  names  does  hot 
constitute  an  official  endorsement  or  approval  of  the  use 
thereof. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 


Army  Research  Laboratory 

Adelphi,  MD  20783-1197 

ARL-MR-494 _ June  2001 

A  One-Dimensional 
Atmospheric  Boundary  Layer 
Model:  Intermittent  Wind 
Shears  and  Thermal  Stability  at 
Night 


Arnold  Tunick 

Computational  and  Information  Sciences  Directorate 


Approved  for  public  release;  distribution  unlimited. 


Abstract 


A  one-dimensional,  time-dependent  computer  model  of  the 
atmospheric  boundary  layer  was  developed  to  simulate  intermittent 
turbulence  and  the  near-ground  microclimate  under  nighttime  stable 
conditions.  In  this  study,  the  model  produced  several  turbulent 
events  (oscillations)  through  the  nighttime  period  that  varied  in 
number,  frequency,  and  strength  along  the  axes  of  initial  geostrophic 
wind  speed.  These  results  were  found  to  be  in  close  agreement  with 
results  from  several  previous  observational  and  theoretical  studies 
of  this  type.  It  is  suggested,  therefore,  that  the  one-dimensional 
computer  model  is  a  useful  mathematical  representation  of  the 
nighttime  case  that  includes  intermittency. 
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1.  Introduction 


Field  experimenters  have  frequently  observed  turbulent  "bursts"  in  the 
lower  boundary  layer,  suggesting  that  they  are  likely  the  result  of  gravity- 
wave  activity  and  other  wind-shear  or  thermal  instabilities  (e.g.,  Schubert, 
1977;  Lu  et  al,  1983;  Coulter,  1990;  Nappo,  1991;  Weber  and  Kurzeja,  1991; 
and  Blumen  et  al,  1999).  Such  turbulent  events  have  even  been  linked  to 
particular  mesoscale  events  and  synoptic  situations,  like  outflow  from  nearby 
thunderstorms,  as  reported  in  Zhou  Ming-yu  et  al  (1980)  and  Neff  (1980). 
In  contrast,  atmospheric  computer  modelers  (e.g.,  Blackadar,  1979;  McNider 
and  Pielke,  1981;  Lin,  1990;  ReVelle,  1993;  and  McNider  et  al,  1995)  describe 
time-dependent  oscillations  at  night  as  resulting  from  nonlinear  interac¬ 
tions  between  the  ground  surface  and  the  airflow  aloft.  From  the  modeler's 
point  of  view,  several  external  parameters,  in  the  right  combination,  can 
increase  the  chances  for  a  model  to  produce  intermittent  turbulence  during 
the  nighttime  period.  These  parameters  are  typically  surface  roughness, 
initial  soil  moisture  or  heat  capacity,  and  geostrophic  (upper  level)  wind 
speed.  Since  intermittent  "bursts"  of  turbulence  in  the  lower  boundary  layer 
cause  wind  gusts  and  fluctuations  in  temperature,  much  attention  has  been 
directed  toward  nighttime  cases  that  include  intermittency  in  modeling  at¬ 
mospheric  pollution  (Nappo  and  Bach,  1997;  Nappo  and  Johansson,  1998). 

Model  results  for  the  nighttime  case  that  includes  intermittency  have  been 
shown  in  previous  works,  e.g.,  Lin  (1990)  and  ReVelle  (1993),  along  with 
descriptions  of  the  parameter  space  within  which  time-dependent  oscilla¬ 
tions  occur.  From  these  earlier  discussions,  however,  it  is  not  clear  to  what 
extent  variations  initial  geostrophic  wind  speed,  surface  roughness,  and 
soil  moisture  affect  the  onset  and  cycling  of  wind  shear  and  thermal  insta¬ 
bilities  at  the  surface  and  aloft.  In  this  report,  therefore,  the  author  presents 
an  expanded  and  improved  view  of  the  nighttime  intermittently  stable  so¬ 
lution  within  the  continuum  of  several  varying  model  parameters.  The  au¬ 
thor  provides  incremental  time  series  of  potential  temperature  and  wind 
speed  as  well  as  time-height  series  of  the  derived  momentum  and  heat  flux 
for  the  layers  2  <  z  <  150  m . 
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2.  The  Atmospheric  Boundary-Layer  Model 


The  one-dimensional,  time-dependent  model  used  in  this  study  is  based  on 
the  previous  works  of  Pielke  and  Mahrer  (1975),  McNider  and  Pielke  (1981), 
and  Avissar  and  Mahrer  (1988).  It  is  a  first-order,  local-turbulence  closure 
computer  model  formulated  to  calculate  the  transfer  of  momentum,  heat, 
and  moisture  at  the  surface  and  aloft.  The  model  uses  an  implicit  finite 
difference  scheme  to  integrate  the  boundary  layer  and  soil-diffusion  equa¬ 
tions  and  contains  many  detailed  formulations  for  the  surface-energy  bud¬ 
get.  The  model's  surface-layer  turbulence  scaling  is  as  described  by 
Zilitinkevich  (1970)  and  Businger  et  al  (1971).  A  formulation  suggested  by 
Smeda  (1979)  is  used  for  the  time-dependent  calculation  of  the  depth  of  the 
nighttime  stable  layer.  Table  1  summarizes  the  model  parameters  for  this 
study  of  the  lower  nighttime  boundary  layer. 


Table  1.  Nighttime 
boundary-layer  model 
parameters. 


Number  of  vertical  levels 

19 

Latitude,  longitude 

(2, 5, 15,  30,  80, 150, 208, 445,  690, 
943, 1203, 1470, 1743,  2023, 2311, 
2609, 2916,  3250,  3500  m) 

45.0  N,  0.0  E 

Surface  roughness 

0.28  m 

Canopy  height 

2.0  m 

Surface  emissivity 

0.95 

Soil  water  content 

0.1  m3/m3 

Average  soil  density 

1400  kg/m3 

Soil  texture 

28.0,  70.0, 2.0 

(%  sand,  %  clay,  %  organic) 

Geostrophic  wind  speed 

u  =  6.5,  8.0, 9.5  m/s  ;  v0  =  0.0  m/s 
91  g 

Day  of  the  year 

Initial  time  (sunset) 

~18  local  time 

Initial  potential  temperature 

-284.0  K 

in  the  surface  layer 

Time  step 

3s 
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The  model  equations  for  the  time-dependent  calculation  of  the  winds  ( u - 
and  f-components),  potential  temperature  (9),  and  specific  humidity  (</)  over 
flat  earth  can  be  expressed  as 


(1) 


dv 

~dt 


d 


") 


K„ 


dzj 


(2) 


96 


d_f 

dz 


K  96 


and 


(3) 


K=9-(k  <*L) 

dt  “  dz.) 


(4) 


where/denotes  the  Coriolis  parameter,  subscript  g  refers  to  the  geostrophic 
wind,  Km  denotes  the  turbulence  transfer  coefficient  for  momentum,  and  Kh 
=  Kq  denotes  the  transfer  coefficients  for  heat  and  moisture.  In  the  surface 

layer,  Km  and  Kh  are  calculated  as  Km  =  ku'^  and  Kh  =  ku'^  ,  respec¬ 
tively,  where  k  is  Karman's  constant,  z  is  height  above  ground  level  (agl)  (in 
meters),  ut  is  the  friction  velocity  (in  units  of  m/ s),  and  <pm  and  (j)h  are 
nondimensional  lapse-rate  functions  that  account  for  surface-layer  stabili¬ 
ties  other  than  neutral.  A  list  of  symbols  is  provided  in  the  appendix. 

During  the  nighttime,  when  the  surface  layer  is  generally  stable,  that  is,  6t  > 

0,  where  a  =  —  —  is  the  potential  temperature  scaling  constant,  the  model 

<t>i,  * 

derives  the  diffusion  coefficients  above  the  surface  layer  as  suggested  by 
Blackadar (1979) 


Km(z)  =  K„(z)  =  \ 


sl2(\  - 18/?/)°'5,  Ri<  0 

si2  (Ricri,  -  Ri)/Ricrit,  0  <Ri<  Ricril 


0,  Ri  >  Ricri, 


(5) 


where  s  is  the  local  wind  shear.  The  parameter,  Ri,  called  the  Richardson 
number,  is  the  ratio  of  thermal  to  mechanical  (wind  shear)  production  of 

turbulent  energy,  so  that  Ri  -  ^  ^  j du/ dz)  +  ( dv/dz )  j.  Riclit  is  the  limit¬ 
ing  value  of  the  Richardson  number,  where  it  is  often  assumed  that 
Ricri,  =  0.25.  Alternatively,  Ricrit  =  0.115Aza175  (McNider  and  Pielke,  1981; 
Avissar  and  Mahrer,  1988),  where  Az  is  the  model  grid  spacing  in  centime¬ 
ters.  This  expression  results  in  having  values  approaching  0.25  for  finer 
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grid  spacing  and  values  closer  to  0.65  where  the  grid  is  more  coarse  (i.e., 

150  to  250  m).  The  length,  /  (in  meters),  /  =  kz(  1  + - — - )  is  eener- 

^  0.0063  w,//  J  6 

ally  thought  of  as  the  width  of  turbulence  and  is  characterized  by  this  for¬ 
mulation  reported  in  Blackadar  (1979),  for  z  above  the  surface  layer.  As 
mentioned  earlier,  the  formulation  suggested  by  Smeda  (1979)  for  the  depth 
of  the  nighttime  planetary  boundary  layer,  zir  is  calculated  as 


(6) 


Some  previous  results  obtained  with  the  use  of  this  expression  for  the  depth 
of  the  nighttime  planetary  boundary  layer  were  in  good  agreement  in  com¬ 
parison  to  the  observed  data  (Tunick,  2000). 

The  surface-energy  budget  as  described  by  Pielke  (1984)  and  Avissar  and 
Mahrer  (1988)  can  be  written  as  . . 


0  ~  A)R,l  +  RlI  ~  RlT  ~  Pc/''a  “  +  Qs  =  F  ,  (7) 

where  A  is  surface  reflectivity  (albedo);  Rsl,  RLi,  and  RL T  are  the  incoming 
solar,  incoming  long-wave,  and  outgoing  long-wave  radiative  fluxes,  re¬ 
spectively,  u.,  0t,  and  cjt  are  the  surface-layer  turbulence  scaling  parameters 
for  wind  speed,  temperature,  and  moisture,  in  that  order;  p  is  the  air  den¬ 
sity;  cp  is  the  specific  heat  of  air  at  constant  pressure;  Lv  is  the  heat  of  trans¬ 
formation  for  water  vapor;  Qs  is  the  soil  heat  flux;  and  F  is  the  function 
applied  in  solving  iteratively  for  the  surface  temperature.  The  formulation 
for  downward  long-wave  radiation  is  based  on  the  earlier  works  reported 
by  Paltridge  and  Platt  (1976),  which  can  be  expressed  as 

RlI  =  -170.9  +  1.195  a  T*  +  0.3  cld£c  o  T4C  ,  (8) 

where  a  =  5.6697  x  10-8  W  nr2  deg-4,  Tr  is  the  reference  level  (~2  m)  tem¬ 
perature  in  kelvins,  cld  is  the  cloud  amount,  ec,  is  the  emissivity  of  the  cloud 
base,  and  Tc  is  the  temperature  of  the  cloud  base  in  kelvins.  This  expression 
has  worked  quite  well  in  two  previous  studies  (i.e.,  Rachele  and  Tunick, 
1994;  Tunick  et  al,  1994). 

Integration  of  the  partial  differential  equations  (1)  to  (4)  is  achieved  with  a 
generalized  form  of  the  Crank-Nicholson  implicit  finite  difference  scheme 
(Paegle  et  al,  1976;  also  see  Pielke,  1984,  sect.  10.1.2,  and  Ahlberg  et  al,  1967). 
Boundary  conditions  applied  at  the  top  of  the  model  are  such  that  the  vari¬ 
ables  retain  their  initial  values. 


3.  Model  Results 


As  reported  in  Businger  (1973),  McNider  et  al  (1995),  and  Mahrt  (1998),  the 
following  events  frequently  characterize  the  nighttime  stable  boundary  layer: 


(a)  Strong  surface  radiative  cooling  and  the  formation  of  a  stably  strati¬ 
fied  layer  near  the  ground. 

(b)  Suppressed  turbulence  as  the  ground  cools  and  the  development  of  a 
laminar  layer. 

(c)  Heat  and  momentum  transfer  of  from  higher  layers  is  blocked  and 
the  Richardson  number  exceeds  its  critical  value  at  some  height  above 
the  ground. 

(d)  Surface  wind  speeds  decrease  (calm).  The  nocturnal  low-level  jet  in¬ 
tensifies. 

(e)  Above  the  laminar  layer,  momentum  continues  to  be  transferred  down¬ 
ward. 

(f)  Above  the  laminar  layer,  little  heat  is  exchanged  since  the  thermal 
gradient  at  this  level  is  close  to  neutral. 

(g)  If  there  is  very  strong  thermal  stability  at  the  surface,  the  airflow  aloft 
may,  in  effect,  become  disconnected  or  decoupled  from  it. 

(h)  Winds  above  the  laminar  layer  will  continue  to  accelerate  in  the  ab¬ 
sence  of  surface  friction. 

(i)  Temperatures  at  the  surface  will  continue  to  decrease. 

(j)  Wind  shears  above  the  laminar  layer  will  build  up  in  the  absence  of 

offsetting  heat  flux. 

(k)  The  Richardson  number  begins  to  decrease  below  its  critical  value. 
The  laminar  surface  layer  will  be  eroded  from  above  by  wind  shears. 

(l)  Downbursts  of  turbulence  cause  wind  gusts  and  oscillations  in  tem¬ 

perature  (and  moisture).  Vertical  mixing  depletes  strong  wind  shears. 

(m)  The  series  of  events  (a)  through  (1)  repeats  itself  through  the  night¬ 
time  period.  This  is  the  intermittently  stable  nighttime  case. 

(n)  If,  however,  winds  and  wind  shears  are  initially  strong  through  the 
lower  nighttime  boundary  layer,  a  more  turbulently  mixed,  weakly 
stratified  layer  is  maintained. 

(o)  Windier  and  warmer  conditions  at  the  surface  will  prevail.  This  is  the 
coupled  nighttime  boundary-layer  solution. 


Figure  1  shows  the  coupled,  decoupled,  and  intermittently  stable  nighttime 
results  for  (a)  potential  temperature  and  (b)  modeled  wind  speed  in  the 
surface  layer  at  5  m  agl.  The  dashed,  circled,  and  solid  lines  are  for  levels  of 
stability  corresponding  to  three  different  initial  values  of  geostrophic  wind 
speed,  i.e.,  6.5,  8.0,  and  9.5  m/ s,  respectively.  Where  turbulence  is  main¬ 
tained  by  strong  wind  shears  (dashed  line),  temperatures  and  wind  speeds 
in  the  surface  layer  remain  high  and  highly  coupled  to  the  air  aloft.  Where 
turbulence  is  suppressed  by  strong  thermal  stability  (solid  line),  the  model 
produces  lower  surface  temperatures*  and  lower  wind  speeds.  The  inter¬ 
mittently  stable  solution  (circled  line)  clearly  oscillates  between  the  coupled 
(windier,  warmer)  and  decoupled  (calmer,  cooler)  states.  Figure  1  of  the 


potential  temperature 
(K)  at  5  m  agl  in  surface 
layer  (a)  and  wind 


values  of  geostrophic 


line). 


*  Certainly,  the  temperature  results  in  figure  1  are  unrealistically  low.  The  so-called  "run¬ 
away"  cooling  (as  expressed  by  Mahrt,  1998)  is  apparently  a  common  defect  in  numerical 
models  that  calculate  meteorological  gradients  near  the  ground  at  night  with  the  use  of 
Obukhov  similarity  functions  (Monin  and  Obukhov,  1954).  In  this  situation,  not  enough 
heat  is  transferred  from  aloft  to  offset  steady  radiative  and  surface  losses. 
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present  study  looks  very  similar  to  ReVelle's  (1993)  figures  3  and  4,  where, 
over  a  period  of  16  hours,  four  oscillations  in  temperature  and  wind  speed 
were  shown  to  occur  at  approximately  4.0, 7.5, 10.5,  and  13.0  h  past  sunset. 
In  this  example,  three  oscillations  occurred  (over  a  period  of  12  hours)  at 
approximately  3.75, 6.5,  and  10.0  h  past  sunset,  which  is  in  close  agreement 
with  ReVelle's  results.  By  this  example,  the  model  is  also  found  to  be  in 
close  agreement  with  several  other  observational  and  theoretical  results  (re¬ 
ported  in  the  literature)  on  nighttime  downburst  events.  Table  2  gives  a 
comparative  summary  of  the  current  and  previous  results. 


Figure  2  contains  data  from  several  consecutive  model  runs  (i.e.,  fig.  2  con¬ 
tains  fig.  1).  The  incremental  time  series  of  potential  temperature  (a)  and 
wind  speed  (b)  in  the  surface  layer  at  5  m  agl  show  that  several  oscillations 
are  generated  through  the  nighttime  period  that  vary  greatly  in  number, 
frequency,  and  strength  along  the  axes  of  initial  geostrophic  wind  speed. 
Figure  2  shows  that  the  onset  of  such  events,  in  hours  past  sunset,  also  var¬ 
ies  along  this  dimension.  The  numbers  of  oscillations  that  occur  are  shown 
to  increase  from  1  to  8  or  more  in  the  direction  of  increasing  of  geostrophic 
winds.  At  the  same  time,  as  the  total  numbers  of  oscillations  increase,  their 
intensity  or  strength  tends  to  decrease.  Figure  2  also  indicates  where  the 
decoupled  nighttime  solution  occurs  (i.e.,  an  extreme  depression  in  tem¬ 
perature  and  calmer  winds,  as  explained  above).  In  contrast,  the  coupled 
nighttime  boundary-layer  solution  is  indicated  where  temperatures  decrease 
less  severely.  An  unexpected  result  for  this  case  is  that  the  model  produces 
oscillations  in  two  cycles.  This  result  is  interesting  because  it  suggests  that  a 
layer  with  sufficiently  strong  wind  speeds  and  wind  shears  can  become 
intermitently  turbulent  at  even  higher  windspeeds  under  the  right  set  of 
conditions. 


Table  2.  Summary  of  previous  observational  and  theoretical  results  on  nighttime  "bursting"  events. 


Source  a 

Type 

(0/T)b 

Surface 

roughness 

(m) 

Geostrophic  wind 
speed 
(m/s) 

Average  number 
of  events  per  night 

Thorpe  and 

O 

not  indicated 

<10.0 

3  hr1  per  event 

Guymer  (1977) 
Schubert  (1977) 

o 

~1.5 

not  indicated 

1-2  within  -5  h 

Blackadar  (1979) 

T 

1.0 

8.0 

-  5-6 

Coulter  (1990) 

O 

0.1-1.0 

not  indicated 

4-5 

Lin  (1991) 

T 

1.0 

8.0 

2-3 

Nappo  (1991) 

O 

0.1-1.0 

not  indicated 

13-18 

Revelle  (1993) 

T 

0.1-1.0 

1.0-3.0 

3-8 

Present  study  c 

T 

0.28  (0.1-1.0) 

8.0  (3.0-11.0) 

3  (1-8) 

aAdapted  from  ReVelle  (1993),  table  2,  p.  1177. 
bO-observations;  T-theory. 

cNumbers  inside  the  parentheses  are  implied  from  the  model  results  shown  in  figures  2,  5,  and  6. 
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Figure  2.  Surface  plot  of 
potential  temperature 
(K)  at  5  m  agl  in  surface 
layer  (a)  and  wind 
speed  (m/s)  at  5  m  agl 
in  surface  layer  (b)  as  a 
function  of  time  past 
sunset  and  value  of 
initial  geostrophic  wind 
speed. 


Figures  3  and  4  show  times  series  of  potential  temperature  data  and  wind 
speed  data,  respectively,  at  15, 30, 80,  and  150  m  agl.  Away  from  the  surface 
(z  >  30  m),  temperatures  and  wind  speeds  are  shown  to  change  very  little 
through  the  nighttime  period  when  geostrophic  wind  speeds  are  low  and 
heat  transfer  is  limited.  In  contrast,  when  initial  geostrophic  wind  speeds 
are  higher  and  turbulent  transfers  of  heat  and  momentum  are  maintained, 
potential  temperatures  are  shown  to  decrease  steadily  through  the  night¬ 
time  period.  Figures  3  and  4  also  show  that  oscillations  are  affected  through 
the  entire  layer  (z  <  150  m),  although  additional  evidence  regarding  the 
onset  of  the  instabilities  would  be  desirable. 

Figure  5  is  similar  to  figure  2  except  that  it  shows  the  model  results  as  they 
are  affected  by  changes  in  surface  roughness,  that  is,  z0  =  0.01  m  (fig.  5(a) 
and  5(c))  and  Zq  =  1.0  m  (fig.  5(b)  and  5(d)).  (In  calculating  the  data  in  figure 
2,  z0  =  0.28  m.)  These  results  show  that  the  surface  layer  tends  toward  warmer 
temperatures  (approximately  8  to  10  K)  and  calmer  wind  speeds  (approxi¬ 
mately  2  to  3  m/s)  as  surface  roughness  increases  from  0.01  to  1.0  m.  This 
is  likely  due  to  increased  surface  stress  (Munn,  1966).  The  results  also  show 
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Figure  3.  Surface  plot 
of  potential  temper¬ 
ature  (K)  at  (a)  15, 

(b)  30,  (c)  80,  and 
(d)  150  m  agl  in  lower 
boundary  layer  as  a 
function  of  time  past 
sunset  and  value  of 
initial  geostrophic 
wind  speed. 


Figure  5.  Surface  plot  of  potential  temperature  (K)  at  5  m  agl  in  surface  layer  and  wind  speed  (m/ s)  at  5 
m  agl  in  the  surface  layer  as  they  are  affected  by  changes  in  surface  roughness  (i.e.,  z0  =  0.01  m  (a)  and  (c) 
and  z0  =  1.0  m  (b)  and  (d). 

that  turbulence  is  generated  about  1.0  to  2.0  hours  earlier  and  at  1.0  to  2.0 
m/s  lower  geostrophic  wind  speeds  as  surface  roughness  increases.  In  other 
words,  oscillations  are  more  likely  to  occur  over  more  rough  (aerodynami- 
cally  turbulent)  surfaces  with  lower  geostrophic  wind  speeds.  The  distri¬ 
bution  and  frequency  of  the  oscillations  also  varies  with  changes  in  the 
roughness  parameter.  The  numbers  presented  in  parentheses  in  table  2  re¬ 
flect  these  additional  results. 

Figure  6  is  also  similar  to  figure  2  except  that  it  shows  the  model  results  as 
they  are  affected  by  changes  in  initial  soil  moisture,  that  is,  ms  =  0.05  m3/ m3 
(fig.  6(a)  and  (c))  and  ms  =  0.20  m3/ m3  (fig.  6(b)  and  (d)).  (In  calculating  the 
data  in  figure  2,  ms  =  0.10  m3/m3.)  These  model  results  show  that  the  sur¬ 
face  layer  tends  toward  warmer  temperatures  (approximately  2  to  3  K)  as 
soil  moisture  increases.  This  is  likely  due  to  decreasing  radiative  losses  and 
other  surface-energy  budget  considerations  (Munn,  1966).  The  results  also 
show  that  the  intensity  of  the  oscillations  increases  for  higher  initial  values 
of  soil  moisture,  which  is  also  an  effect  related  to  the  heat  capacity  of  the 
soil  and  surface  cooling  (see  Lin,  1990). 
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Figure  6.  Surface  plot  of  potential  temperature  (K)  at  5  m  agl  in  the  surface  layer  (a)  and  (b)  and  wind 
speed  (m/ s)  at  5  m  agl  in  surface  layer  (c)  and  (d)  as  they  are  affected  by  changes  in  initial  soil  moisture 
(i.e.,  ms  =  0.05  m3/m3  (a)  and  (c)  and  ms  =  0.20  m3/m3  (b)  and  (d). 


Figure  7(a),  (b),  and  (c)  shows  the  calculated  momentum  flux,  —u'w'  — 

in  units  of  m2/s2,  and  heat  flux,  -VW  =  Kh(de/dz ) 

in  figure  7(d),  (e),  and  (f)  in  units  of  ms  3  K,  in  the  lower  boundary  layer  for 
levels  of  stability  corresponding  to  three  different  initial  values  of  geostrophic 
wind  speed,  that  is,  6.5,  8.0,  and  9.5  m/ s.  In  the  very  stable  layer  case  (fig. 
7(a)  and  7(d)),  momentum  increases  through  the  lower  boundary  layer  as 
the  nocturnal  low-level  jet  forms  aloft.  The  model  results  show  that  the 
momentum  flux  continues  downward  through  the  lower  boundary  layer, 
even  as  the  surface  layer  becomes  increasingly  more  stable  (and  eventually 
laminar).  This  may  explain  why,  in  figure  1,  the  surface-layer  wind  speeds 
are  around  2.5  m/ s,  rather  than  calm,  through  the  nighttime  period.  In 
contrast,  the  heat  flux  through  the  laminar  layer  decreases  as  surface-layer 
temperatures  (not  shown)  approach  a  kind  of  radiative  equilibrium  with 
the  ground  (Mahrt,  1998).  Above  the  laminar  layer  (z  >  30  m),  the  heat  flux 
is  nearly  absent. 
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Figure  7.  Contour  plot  of  calculated  momentum  flux  (a),  (b),  and  (c)  and  heat  flux  (d),  (e),  and  (f)  through 
nighttime  period  for  values  of  geostrophic  wind  speed  of  6.5  (a)  and  (d),  8.0  (b)  and  (e),  and  9.5  m/ s  (c) 
and  (f). 

In  the  intermittently  stable  case  (fig.  7(b)  and  7(e)),  a  series  of  sporadic  break¬ 
downs  are  shown  as  fluctuations  of  momentum  and  heat  flux  in  the  layer, 
2  <  z  <  150  m.  There  appear  to  be  about  5  or  6  such  events  through  the 
nighttime  period.  For  this  example,  the  surface  roughness  and  initial  geo¬ 
strophic  wind  speed  were  0.28  m  and  8.0  m/ s,  respectively  (same  as  fig.  1). 
Further  analysis  of  these  results  would  be  desirable  to  provide  additional 
information  regarding  the  onset  wind  shear  and  thermal  instabilities  aloft. 

In  the  weakly  stable  case  (fig.  (c)  and  (f)),  both  the  heat  and  momentum  flux 
remain  strong  through  the  nighttime  period.  The  increased  wind  speeds 
(and  wind  shears)  in  this  example  do  not  permit  the  formation  of  a  sepa¬ 
rate,  highly  stratified  layer  at  the  surface.  Instead,  the  model  derives  a  steady 
and  much  less  severe  decrease  in  temperatures  (not  shown)  through  the 
entire  layer,  that  is,  2  <  z  <  150  m.  This  is  the  coupled  nighttime  boundary- 
layer  result. 
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4.  Conclusion 


In  this  study  a  computer  atmospheric  model  is  developed  to  simulate  the 
lower  nighttime  stable-boundary  layer.  An  overview  is  given  for  the  devel¬ 
opment  of  wind  shear  and  thermal  stability  in  the  lower  boundary  layer 
under  weakly  stable,  highly  stable,  and  intermittently  stable  conditions. 
Several  graphs  show  incremental  time  series  of  the  model  data.  They  pro¬ 
vide  an  improved  and  expanded  view  of  the  nighttime  wind  and  tempera¬ 
ture  oscillations  and  of  the  parameter  space  within  which  such  time-depen¬ 
dent  oscillations  occur.  The  model  used  in  this  study  is  shown  to  be  a  useful 
representation  of  the  nighttime  case  that  includes  intermittent  turbulence. 
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Appendix.  Symbols 


A  surface  reflectivity  (albedo) 

cld  cloud  amount  (in  tenths) 

c  specific  heat  of  air  at  constant  pressure 
f  Coriolis  parameter 

F  function  applied  in  solving  for  the  surface  temperature  through  the 
surface-energy  budget 
g  acceleration  due  to  gravity 

k  Karman's  constant  (~0.4) 

Kh  eddy  transfer  coefficient  for  heat 

Km  eddy  transfer  coefficient  for  momentum 

K  eddy  transfer  coefficient  for  moisture 

/  mixing  length 

Lv  heat  of  transformation  for  water  vapor 
ms  soil  moisture  or  soil  water  content 
q  specific  humidity 

qt  surface-layer  turbulence  scaling  parameter  for  moisture 

Qs  soil  heat  flux 

Ri  ratio  of  thermal  to  mechanical  (wind  shear)  production  turbulent 
energy  called  the  Richardson  number 
Ricrjt  limiting  value  of  the  Richardson  number 
Rsl  incoming  solar  radiative  flux 

RLI  incoming  long-wave  radiative  flux 

RLI  outgoing  long-wave  radiative  flux 

s  local  wind  shear, 

t  time  in  hours  past  sunset 

Tc  temperature  of  the  cloud  base  (in  kelvins) 

T  reference  level  (~2  m)  temperature  (in  kelvins) 

u  east-west  component  of  horizontal  wind  speed 

Ug  east-west  component  of  geostrophic  wind  speed 
«,  surface  friction  velocity 

u'  eddy  fluctuating  component  of  the  horizontal  wind  speed 
v  north-south  component  of  horizontal  wind  speed 
v  north-south  component  of  geostrophic  wind  speed 
vf  eddy  fluctuating  component  of  the  vertical  wind  speed 
£  emissivity  of  the  cloud  base 

<ph  nondimensional  temperature  lapse  rate 

(j)m  nondimensional  wind  shear 

6  potential  temperature 

8  potential  temperature  scaling  constant 

p  air  density 
a  Stefan-Boltzmann  constant 

any  overbar  denotes  the  mean 
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